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I. INTRODUCTION 

Strongly correlated one-dimensional (ID) quantum liq- 
uids and gases are nowadays available in a large number 
of different laboratory systems ranging from single-wall 
carbon nanotubes 1 to semiconductor nanowires 2 ", con- 
ducting molecules^, and trapped atomic gases^Si 7 .. Chi- 
ral Luttinger liquids at fractional quantum Hall edges^ 
also provide a beautiful example of ID conducting quan- 
tum liquids and have been the subject of intense experi- 
mental 9 and theoretical efforts^ -. 

There are two fundamental key features that are com- 
mon to all these ID systems, (i) Independently of statis- 
tics, their effective low-energy description is based on a 
harmonic theory of long- wavelength fluctuations^ due to 
the interplay between topology and interactions, (ii) In 
the most interesting and exciting experimental situations 
the translational invariance of the liquid can be broken 
due to the presence of inhomogeneous external fields of 
different types, such as magnetic traps in the case of ul- 
tracold atomic gases 1 ^ or Hall bar constrictions in the 
case of fractional quantum Hall edges^. These strong 
perturbations induce the appearance of new length scales 
causing novel physical behaviors relative to the corre- 
sponding unperturbed, translationally invariant model 
system. 



A powerful theoretical tool to study the inter- 
play between interactions and inhomogeneous external 
fields of arbitrary shape is density-functional theory 
(DFT) 13 i 14 i 15 i 16 , which is based on the Hohenberg-Kohn 
theorem 1 ^ and on the Kohn-Sham mapping to an aux- 
iliary noninteracting system 1 —. Many-body effects en- 
ter DFT via the exchange-correlation (xc) functional, 
which is often treated by the local-density approxima- 
tion (LDA) 13 i 14 i 15 i 16 i 17 i 18 . The essence of LDA is to lo- 
cally approximate the xc energy of the inhomogeneous 
system under study with that of an interacting homoge- 
neous reference fluid, whose correlations are transferred 
by the LDA to the inhomogeneous system. For exam- 
ple, for inhomogeneous 2D and 3D electronic systems the 
underlying reference fluid is normally the homogeneous 
electron liquid (EL) 1 ^ 1 ^, whose xc energy is known to a 
high degree of numerical precision thanks to the quan- 
tum Monte Carlo (QMC) technique 20 . However, these 
ELs are believed to be normal Fermi liquids 1 ^* 1 ^ over a 
very broad range of densities, whereas their ID analogue 
is described by the Luttinger-liquid model-i. Thus inho- 
mogeneous II? fermionic systems appear as an interesting 
example in which it is appropriate to change the refer- 
ence system to one that possesses ground-state Luttinger- 
liquid rather than Fermi-liquid-type correlations22iSi24. 

Several other examples have been discussed in the lit- 
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erature in which either the LDA reference system is not 
an EL or the auxiliary system of the Kohn-Sham map- 
ping is not an assembly of noninteracting particles. Kohn 
and co-workers have introduced the concept of the "edge 
electron gas"— to study electronic edge regions where 
the single-particle wavefunctions evolve from oscillatory 
to evanescent. In the presence of broken gauge symme- 
try, such as in the superconducting state, an appropriate 
reference system is the EL in the presence of an external 
pairing field inducing superconducting correlations^^ 7 .. 
A related approach has also been proposed for Bosc- 
Einstcin-condcnsed systems^. Similar in spirit to the 
above mentioned work on DFT for the Hubbard model, 
is DFT for the Heisenberg model, in which the reference 
system is a lattice of spins on equivalent sites^. Finally, 
we mention spectral-DFT— , in which the Kohn-Sham 
noninteracting system is replaced by a suitably chosen 
interacting system not handled via the usual LDA. Such 
departures from the standard EL paradigm substantially 
expand the range of usefulness of DFT in condensed- 
matter physics, but also demand the construction and 
investigation of new classes of functionals. The present 
work is concerned with the testing of LDA-type density 
functionals for strongly correlated II? ultracold Fermi 
gases confined inside an optical lattice. 

From the experimental point of view such systems, 
which arc highly tunable and ideally clean, are attracting 
a great deal of interdisciplinary interest because they al- 
low to realize strongly interacting many-body systems 
through the manipulation of relevant degrees of free- 
dom other than the bare atom-atom interaction, such 
as the well depth of the optical lattice that allows to 
tune the relative strength of hopping to on-site repul- 
sion/ attraction^!. In particular the study of these sys- 
tems may help us understand a number of phenomena 
that have been predicted in solid-state and condensed- 
matter physics. Several effects, known in these subficlds 
of physics for decades, have already been observed and 
quantitatively analyzed in ultracold atomic gases trapped 
in optical lattices. Two beautiful examples are the Bloch 
oscillations under an applied force in a ID optical lat- 
ticed and the superfluid-to-Mott insulator transition of 
a Bose-Einstein condensate in a 31? optical lattic e - 33 ! 34 . 
Typical ID quantum phenomena have also been observed 
in both Bose and Fermi gases. For instance, in the 
work of Paredes et al~ and of Kinoshita et al& a 87 Rb 
gas has been used to realize experimentally a Tonks- 
Girardeau system^. The more recent preparation of two- 



component 40 K Fermi gases in a quasi-l-D geometry^ pro- 
vides a unique possibility to experimentally study phe- 
nomena that were predicted a long time ago for electrons 
in a ID solid-state environment, such as spin-charge sep- 
aration in Luttinger liquids 1 ^! and charge-density waves 
in Luther-Emery liquids 3 ^. 

In this work we focus our attention on a particular 
ID lattice system of ultracold atoms: a two-component 
Fermi gas with repulsive intercomponent interactions in 
the presence of static external potentials that break the 
lattice translational invariance. Theoretical studies of 
this model have been carried out both by numerical tech- 
niquea 37 i 38 i 39 and by LDA-based calculationa 40 i 41 i 42 that 
will be discussed in detail later in this work. Building 
upon the earlier ideas described in Refs. I22I23L we here 
employ a lattice DFT scheme in which the xc potential is 
determined exactly at the LDA level through the Bcthc- 
Ansatz solution of the homogeneous ID Hubbard model. 
The results are tested against accurate QMC simulation 
data over a broad range of values for the Hubbard on- 
site interaction, the number of atoms and lattice sites, 
and different types of external potential. 

The contents of the paper are briefly described as fol- 
lows. In Sect. |H] we introduce the lattice Hamiltonian 
that we use to describe the system of physical interest. 
For the benefit of readers who are not familiar with the 
Bethe-Ansatz and the Luttinger liquid, we also briefly 
summarize the properties of the model and its solution in 
the absence of external potentials. In Scct. lIIII we present 
the self-consistent lattice DFT scheme that we use to deal 
with the inhomogeneous system under confinement and 
explain in detail the Bethe-Ansatz LDA that we employ 
for the xc potential. In Sect. II Vl we report and discuss our 
main theoretical results in comparison with QMC simu- 
lation data. Finally, in Sect. we summarize our main 
conclusions. An Appendix contains the formal derivation 
of the lattice Kohn-Sham equations. 



II. THE FERMIONIC HUBBARD MODEL 

We consider a two-component interacting Fermi gas 
with Nf particles which are constrained to move under 
confinement inside a ID lattice with unit lattice constant 
and N s lattice points labeled by the discrete coordinate 
Zi = i,i G [1, N s ], This system is described by the follow- 
ing single-band Hubbard Hamiltonian^M 3 -, 



i,j cr i i 

1 



where Uj = t > if i and j are nearest-neighbor sites and degree of freedom (hypcrfine-state label) . The field op- 
zero otherwise, and a =|, { represents a pseudospin-1/2 
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erator c\(zi) (c a (zi)) creates (destroys) a fermion with 
pseudospin a at position z i} n a (zi) — c^ r (zi)c <T (Zi) is 
the pseudospin-resolved site occupation operator nor- 
malized to the number of particles with pseudospin a, 

N(T = {J2i "<r( z i))> an d fl(zi) = J2a^( Z i) is tne total 

site occupation operator with (X^"-( z 0) = -^f- Finally 
Vext(zi) is an external static potential associated with the 
confinement. The Hubbard Hamiltonian without con- 
finement has, of course, been widely used in studies of 
strongly correlated electrons, where the site index refers 
to ion positions. It also applies to an electron liquid in a 
quantum wire under the effect of a spatially modulated 
electric potential^, where the site index refers to min- 
ima of the superlattice structure. Parabolic confinement 
is easily achieved in nanowire quantum dots 4 ^, and most 
of the results reported below also apply to that case. 

The physical processes associated with each term in 
Eq. are clear. T describes kinetic processes of 
atom hopping with site-to-site tunneling amplitude t. 
7iint is the interspecies on-site Hubbard interaction with 
strength U. The intraspecies scattering can be assumed 
to be negligible because atoms with the same pseudospin 
are kept apart by the Pauli principle and can thus be 



taken, to a very good approximation, as noninteracting 
in an ultracold gaseous state. 7i C xt gives the coupling 
of the atoms to a static external potential. In this work 
we will restrict our attention to repulsive interactions, 
i.e. U > 0, in symmetric systems with equal numbers of 
atoms of each pseudospin species (iV-j = iVj = Nf /2). At- 
tractive interactions have been discussed in Ref. |4l| and 
asymmetric Fermi gases in Ref. \a$L 

In the absence of a longitudinal external field (V cxt = 
0), H reduces to the Hamiltonian of a ID homogeneous 
Hubbard model (HHM) that has been solved exactly 
more than 30 years ago by Lieb and Wii 4 ^. At zero 
temperature and for N-\ — N± the properties of the 
ID HHM in the thermodynamic limit (Nf,N s — > oo) 
are determined by two parameters only, the filling fac- 
tor n — Nf/N s and the dimensionless coupling constant 
u = U/t. 

According to Lieb and Wii 4 ^, the ground state (GS) 
of the repulsive ID HHM in the thermodynamic limit is 
described by two continuous distribution functions p(x) 
and cr(y) which satisfy the Bethe-Ansatz (BA) coupled 
integral equations, 
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The parameter Q is determined by the normalization con- 
dition J^q p(x)dx — n, while <r(y) is normalized accord- 
ing to J^/ 00 cr(y)dy = n/2. The GS energy of the system 
(per site) is given by 



3 (n <l,u) = -2t 



+Q 



p(x) cosxdx . 



(4) 



For repulsive interactions the ID HHM describes a Lut- 
tinger liquid 48 if n ^ 1 or 2. At half-filling, i.e. n = 1, the 
GS is a Mott insulator for every u ^ 0, while for n = 2 
it is a band insulator. The two metallic GS branches for 
n < 1 and n > 1 are connected by particle-hole symme- 
try, £ G s(^ > 1, u ) = (n — 1)U + £ G s(2 — n, it). Note that 
the presence of a Mott-insulating GS at half-filling is sig- 
nalled by a cusp in the GS energy at n = 1, induced by 
the linear term in this relation (see Fig.^J. Correspond- 
ingly the charge excitation spectrum possesses a gap. 

The GS energy is analytically known for u = 
(noninteracting fermions) as s GS (n < l,u — 0) = 
—At sin (mr/2)/7T, for u = +oo as £ G s(^ < 1> u — * +oo) = 



— 2t sin (nn) /it, and for every positive value of u at half- 
filling— as 



3 (ri = l,u) = -At 



Jo(x)Ji(x) 
x[l + exp (xu/2)] 



dx . 



(5) 



In Fig.^ we show e G s as a function of n for various values 
of u. The inset in Fig. ^ shows the cusp at n = 1. 

In the next Section we show how the GS properties 
of the inhomogeneous system described by Eq. (JTJ can 
be calculated very accurately in a DFT scheme based on 
Eqs. 



III. LATTICE DENSITY-FUNCTIONAL 
THEORY 

A powerful tool to calculate the GS properties of a 
Hamiltonian such as given in Eq. Q is a lattice version 
of DFT, the so-called site-occupation functional theory 
(SOFT). This was introduced in pioneering papers by 
Gunnarsson and Schonhammer— to study the band-gap 
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FIG. 1: Ground-state energy of the repulsive ID homoge- 
neous Hubbard model (in units of hopping parameter t) as a 
function of the filling factor n for various values of the cou- 
pling parameter u. The cusp at n = 1 signals the Mott- 
insulating phase. The inset shows an enlargement of the re- 
gion 0.95 < n < 1.05. 



problem in the context of ab initio theories of fundamen- 
tal energy gaps in semiconductors and insulators^. We 
summarize the key theoretical aspects of SOFT in the 
Appendix ^ in order for the present paper to be self- 
contained. 

Within SOFT the exact GS site occupation n GS (zi) — 
(GS|n(zi)|GS) can be obtained by solving self- 
consistently the lattice Kohn-Sham (KS) equations 



£gs [^gs] 



''xc (^i)^GS (Zj) 



J2Unl s ( Zl )/4 + £ x 



(8) 



Equations JBJ-lJSJ provide a formally exact scheme to 
calculate n GS (zi) and £ GS , but £ xc needs to be approxi- 
mated. The LDA has been shown to provide an excellent 
account of the GS properties of a large variety of inhomo- 
geneous systemsi2ii2ii£ii£iiLi&. In this work we employ 
the following BA-based LDA (BALDA) functional 

<c ALDA Ks](^) = v^ m (n, u)\ n ^ nog[zi) , (9) 

where, in analogy with ab initio DFT, the xc potential 



.horn 



(n, u) of the ID HHM is defined by 



(n, u) = — 
on 



U 



e GS (n, u) - e GS (n, 0) - ^n 2 



(10) 



Thus, within the LDA scheme proposed in Eqs. (J5J) 
and l|l(J|) . the only necessary input is the xc potential 
of the ID HHM. which is known from its BA solution. 



A. The exchange-correlation potential of the ID 
HHM 



In what follows we propose two alternative ways to 
calculate the xc potential of the HHM. 



+ v KS [n GS ](zi)8ij](p a (zj) = E a ip a (z.i) (6) 
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BALDA/LSOC 



with v KS [n GS ](zi) = Un GS (z l )/2 + v xc (zi) + V ext (zi), to- 
gether with the closure 

n GS (zi) = ^2 r " \Va{zi)\ ■ (7) 

a,occ. 

Here the sum runs over the occupied orbitals and the 
degeneracy factors r Q satisfy the sum rule J2 a ^ a = ^{ ■ 
The first term in the effective Kohn-Sham potential u KS is 
the Hartree mean-field contribution, while v xc [n GS ](zi) = 
S£ xc [n] / Sn(zi)\ GS is the xc potential defined by the deriva- 
tive of the xc energy £ X cM evaluated at the GS site oc- 
cupation [see Eq. (|A9|) ] . Notice that exchange interac- 
tions between parallel-pseudospin atoms have been effec- 
tively eliminated in the Hubbard model Q by restrict- 
ing the model to one orbital per site. Hence parallel- 
pseudospin interactions are not treated dynamically in 
solving the Hamiltonian, but are accounted for implic- 
itly via a restriction of the Hilbert space. To stress 
the analogy of the present work with ab initio applica- 
tions of standard DFT, we nevertheless continue to call 
£ xc [n] and w xc [«gs] the exchange-correlation energy and 
the exchange-correlation potential, but it is understood 
that the exchange contribution to these quantities is ex- 
actly zero. The total GS energy of the system is given 



A semi-analytical scheme, in which the calculation of 
i>x° m (n, it) is carried out with an accurate parametriza- 
tion formula for e GS (n, u), has been proposed by Lima 
et al. (LSOC) 23 . This is very similar in spirit to what 
is done in the EL-based LDA calculations on 3D and 
2D electronic systems*, where the only input is the xc 
energy of the EL for which accurate parametrizations 
are available. Results for n GS (zi) that are obtained with 
v x° m ( n t u ) determined according to this semi-analytical 
route will be labeled by the acronym BALDA/LSOC. 

2. BALDA/FN 

A very appealing feature of Eqs. © and IjlUI) from the 
formal DFT viewpoint, is that one can go a step further 
than the usual parametrized LDA and establish a fully 
numerical improvement over the BALDA/LSOC scheme. 
This procedure does not rely on any approximation for 
v x ° m (n,u) and can easily be set up by observing that 
v x ° m (n,u) satisfies the exact BA equation 

r+Q 

v^° m (n<l,u) = -At / [d n p(x)} cos x dx 
Jo 

- At p(Q)(d n Q) cos Q + 5v KH (11) 
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FIG. 2: The exchange-correlation potential of the repulsive 
ID homogeneous Hubbard model (in units of t) as a function 
of n for various values of u. The inset shows the discontinuity 
A xc (it) (in units of t) as a function of u. 



and the symmetry v^° m (n > l,u) — — v^° m (2 — n,u). In 
Eq. Ijllfl we have Svkh = 2t cos (nir/2) — Un/2. Equa- 
tion Ijlll) must be supplemented by a set of exact BA 
equations for d n p, d n a, and d n Q, which can be derived 
from Eqs. 10- © upon differentiating with respect to 
Illustrative numerical results for v^° m (n,u) are reported 
in Fig. |2 Results for n GS (zi) that are obtained with 
!)J c om (ii,ii) determined according to this fully numerical 
route will be labeled by the acronym BALDA/FN. 

A second appealing feature of the local scheme in 
Eqs. © and ltTU|) is that the Mott cusp in s GS (n, u) is 
responsible for an intrinsic discontinuity A xc (u) in v^° nl 
at n = 1, 

A xc (u) = lim v^° m (n,u)- lim v l ™ m (n,u) 
n— »1+ n— *1~ 

= U-2 lim d n £ GS (n < l,u) (12) 

n— *1 

(see the inset in Fig.©. As a consequence and contrary to 
the EL-based LDA, the xc potential in Eq. (JDJ possesses 
a discontinuity in its derivativeiS^LS. 

Correlation-induced discontinuities, here related to 
Mott-transition physics, also appear in the xc potential 
of a 2D EL in the fractional quantum Hall regime, where 
correlation-induced gaps at fractional filling factors are 
associated with the formation of an incompressible liquid 
(see e.g. Fig. 10.28 in Refii&). These physical discontinu- 
ities make it difficult to obtain converging self-consistent 
solutions of the KS equations whenever the local density 
reaches a value associated with a gap in the homoge- 
neous reference fluidS. In our case this happens only 
when n GS (zi) — > 1 at some position Zj. In the context 
of DFT calculations of the edge structure of fractional 
quantum Hall liquids, Ferconi et alm& have handled this 
convergence problem by going to a very small finite tem- 
perature. This allows fractional occupation of the single- 
particle KS levels. 



3. TLDA 

A conceptually simpler route to solve such convergence 
problems, which we have examined in this work, is to re- 
sort to an LDA also for the noninteracting kinetic energy 
functional %[n] (see the Appendix lATl . This is similar to 
the Thomas-Fermi (TF) approximation, which involves 
an LDA for T s [n] but ignores the exchange and correla- 
tion energy. This in itself is not a fully reliable approach 
for the present work, which is directed at strongly cor- 
related regimes. However, we would like to exploit one 
favorable aspect of the TF approach, i. e. the replacement 
of self-consistent solutions by a direct minimization of to- 
tal energy functionals. To achieve this we combine a TF- 
like LDA for 7^ [n] with an LDA for £ xc [n] . This amounts 
to approximating all nontrivial terms in the total energy 
functional by the LDA, and for this reason we refer to this 
approach as the total-energy LDA (TLDA) 42 . At vari- 
ance from Ref. 0> where TLDA was used in conjunction 
with the LSOC parametrization of the xc energy, we here 
employ its fully numerical counterpart 4 ^, in which all re- 
sults pertaining to the homogeneous reference system are 
obtained from the Lieb-Wu equations. Both formulations 
of TLDA circumvent a self-consistent solution of the KS 
equations, at the expense of a less accurate account of 
the kinetic energy. The TLDA approximation consists in 
writing 

£o T s LD >os] = E £ -MU«cs(z,) 
i 

+ },Vex.t(zi)n GS (zi) . (13) 

i 

Within this approximation the variational equation (|A2|) 
can be written as 

9 n s G s(n,0)\ n ^ nG3 ^ + v KS [n GS ](zi) = constant, (14) 

where the constant is fixed by normalization. Results for 
n GS (zi) that are obtained by this route will be labeled by 
the acronym TLDA. 

IV. NUMERICAL RESULTS 

We now turn to illustrate our main numerical results 
for the density profile of paramagnetic Fermi gases, which 
are summarized in Figs. 1311 II In this work we focus on 
external potentials of the general form2& 

V^(z i ) = ^2Vt{z i -N B /2) t , (15) 
e>i 

where I is an integer and Ve a constant. 

We have solved numerically the self-consistent scheme 
represented by Eqs. ©-© by using both the 
BALDA/LSOC parametrization 23 and the BALDA/FN 
procedure. Results obtained through the TLDA in 
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TABLE I: Ground-state energy (in units of i) of a repulsive 
Fermi gas with Nf = 30 atoms, trapped in a harmonic poten- 
tial of strength V 2 /t = 6 x 10" 3 and in a lattice with N s = 100 
sites. 



u 


BALDA/LSOC 


BALDA/FN 


QMC 


±2 


-17.16 


-16.68 


-16.68 ±0.07 


+4 


-11.33 


-10.84 


-10.74 ±0.08 


+6 


-8.00 


-7.52 


-7.64 ±0.08 


+8 


-5.95 


-5.50 


-5.62 ± 0.11 



Eqs. (|13|) an d (|14f> have also been examined. In par- 
allel, QMC simulations have been performed using a 



zero-temperature projector approac. 



,54,55 



adapted from 



the QMC determinantal algorithm of Scalapino et alm&. 
Within this approach a projector operator exp(— 8H) is 
applied to a trial wave function, which we choose to be 
the exact ground state for u = (we have adopted a 
projector parameter 9 with values up to 45/i for the ac- 
curate comparisons that we present in this work). We 
have used a discrete Hubbard-Stratonovich transforma- 
tion^ in decoupling the fermionic degrees of freedom in 
the interaction term of TL. A detailed description of our 
QMC approach can be found in Refs. 58 59 60. 

In Fig. |31 we show the GS site occupation for a Fermi 
gas with Nf = 30 atoms trapped in a purely harmonic 
potential (i.e. V^2 = 0) of strength Vz/t = 6x 10~ 3 and 
in an optical lattice with N s — 100 sites. The interac- 
tion parameter is increased from u = ±2 to u = ±8. 
The agreement between the BALDA/FN scheme and 
the QMC results is clearly excellent for all values of u. 
The BALDA/LSOC scheme gives similar results and im- 
proves with increasing coupling (it > 6), but overall the 
BALDA/FN results are closer to the QMC data. In the 
rest of the paper we will thus focus on the BALDA/FN 
scheme. 

Computationally, BALDA/FN is slightly more expen- 
sive than BALDA/LSOC. Both typically take a few sec- 
onds on a small PC to generate a single density profile 
like one of those reported in Fig. 01 QMC runs may take 
about 30 hours for a single density profile, but unlike 
LDA also provide access to correlation functions and to 
the momentum distribution. In particular, the calcula- 
tion of the GS energy, which is straightforward within 
the BALDA-DFT scheme, requires a careful extrapola- 
tion procedure within QMC, which is explained in detail 
in Appendix [B] and illustrated in Fig. 0] The extrapo- 
lated QMC GS energies corresponding to the GS density 
profiles shown in Fig. [31 have been reported in Table [I] 
with their estimated statistical error, together with the 
BALDA/FN and BALDA/LSOC results. 

In Fig. we compare TLDA results with BALDA/FN 
and QMC results for the same system parameters as in 
Fig. E| for the cases u = ±2 and u = ±8. TLDA is 
slightly less accurate than BALDA/FN, especially at rel- 
atively small values of the interaction parameter where 
hopping kinetic processes, which are treated at a simple 



LDA level in the TLDA, are still important. For the same 
reason, the regions close to the edge of the trap are those 
where the TLDA is less accurate. 

In Fig. we show the local inverse compressibility, 



L (*i) 



d 2 e c 



i{n,u) 



dn 2 



(16) 



corresponding to the system parameters as in Fig-01 This 
quantity is clearly much more sensitive to the variations 
of the interaction parameter u than the GS site occupa- 
tion. 

In Fig. [7] we test the performance of the BALDA/FN 
scheme upon variations of the strength of the harmonic 
potential V^/t (top panel) and of the particle number 
Nf (bottom panel). Again, the agreement between the 
BALDA/FN scheme and the QMC results is excellent for 
all the values that we have checked. 

We show next some results for the GS density profiles 
of Fermi gases trapped in anharmonic potentials. Fig. [S] 
illustrates two possible situations. In the top panel we 
show an asymmetric external potential which contains 
a main harmonic component with strength Vi/t = 1.6 X 
10~ 2 , a cubic component with strength V^/t = 1.6 x 10 -4 
breaking inversion symmetry, and a quartic component 
with strength Vi/t = 1.92 x 10~ 5 ensuring existence of a 
ground state. The presence of cubic and quartic compo- 
nents represents small deviations from a purely harmonic 
trapping potential that may occur in a real trap in the 
laboratory. In the bottom panel of Fig. |S| we show a 
double-well potential similar to the one that is used to 
model tunncling-coupled lateral semiconductor quantum 
dots: this potential contains a harmonic component with 
strength Vi/t = —4 x 10~ 3 and a quartic component with 
strength V 4 /t = 3 x 10~ 6 . 

In the top panel of Fig.^we show the BALDA/FN pre- 
dictions and the QMC results for the GS site occupation 
of a Fermi gas with Nf — 14 atoms and u = ±8, trapped 
in the asymmetric potential depicted in the top panel of 
Fig. [S] In the same figure TLDA results are also shown. 
The performance of the TLDA scheme at weaker inter- 
actions deteriorates with decreasing particle number. In 
the bottom panel of Fig. [5] we show the GS site occupa- 
tion for a Fermi gas with u = ±2 trapped in the double- 
well potential depicted in the bottom panel of Fig. |SJ In 
Fig. we show the extension of Fig. [5] to the case of at- 
tractive interactions, which can be handled by means of 
the technique described in Ref. l4ll Th e site occupation 
combines features found in Refs. I4ll62l for attractive in- 
teractions in a single parabolic well, i.e. atomic-density 
waves induced by Luther-Emery spin pairing, with fea- 
tures shown in the bottom panel of Fig. [5] for a double 
well and repulsive interactions. Repulsive interactions 
are seen to lead to a higher density at the center of the 
double-well potential, suggesting enhanced tunneling be- 
tween the wells. 

In all these cases the BALDA/FN scheme has been 
found to be extremely accurate, as judged by compar- 
ison to QMC. However, all data shown in Figs. 131 101 
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FIG. 3: Site occupation n GS (zi) as a function of site position Zj for a repulsive Fermi gas with Nf — 30 atoms, trapped in a 
harmonic potential of strength V%jt = 6 X 10~ 3 and in a lattice with N s = 100 sites. The interaction parameter is varied from 
u — +2 in the top left panel to u = +8 in the bottom right panel. BALDA/LSOC and BALDA/FN results are compared with 
the QMC data. 
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correspond to "purely metallic" phases of the interact- 
ing Fermi gas, i.e. phases in which n GS (zi) < 1 every- 
where inside the trap. From earlier work 37,38,40 we know 
that in the trap there can be metallic phases that coex- 
ist with Mott-insulating and/or band-insulating regions, 
i.e. phases in which n GS (zi) is locally locked to 1 or 2. 
What happens with the BALDA/FN scheme if the sys- 
tem parameters are such that n GS (zi) becomes unity or 
very close to unity at some point in space? As we have 
mentioned in Sect. 11111 it is difficult in this case to obtain 
converged self-consistent solutions of the KS equations 
due to the discontinuity at n = 1 in the xc potential of 
the ID HHM. It is still possible to obtain a converged 
self-consistent solution in a very reasonable number of 
iterations if the discontinuity is relatively small. An ex- 
ample is given in the top left panel of Fig. 1111 where we 
show the GS site occupation for a Fermi gas with Nf = 70 
atoms in an optical lattice with N s — 100 sites, subject to 
a purely harmonic potential of strength V^/t = 2.5 xl0~ 3 . 
For u — +2 the size of the xc discontinuity is so small that 
we are able to achieve an accurate self-consistent solution 
of the BALDA/FN equation. Upon further increasing u 
in the same system while keeping unchanged all other 
control parameters, it becomes progressively difficult to 
obtain converged self-consistent solutions of the KS equa- 
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pulsive Fermi gas with N{ = 30 atoms trapped in a harmonic 
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N s = 100 sites. The interaction parameter is u — +2 in 
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BALDA/FN results are compared with QMC data. 
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FIG. 6: BALDA/FN results for the local inverse compressibil- 
ity [7rfv(zi)] -1 (in units of t) as a function of Zi for the same 
system parameters as in Fig. 



tions in a reasonable number of iterations. The easiest 
way to sidestep this problem is to resort to the TLDA 
scheme, Eqs. (|13H and 114(1. which is able to capture the 
main physical features of the above mentioned coexis- 
tence of compressible and incompressible regions^SS. 
In Fig. ^2 we compare TLDA results with QMC re- 



FIG. 7: Top panel: Site occupation n GS (zi) as a function of 
Zi for a repulsive Fermi gas with Nt = 30 atoms, trapped in a 
harmonic potential with strength Vijt = 2 x 10~ 3 ,4 x IO -3 , 
and 6 x 10 -3 (from bottom to top) and in a lattice with N B = 
100 sites. Bottom panel: Site occupation n G s(zi) as a function 
of Zi for a repulsive Fermi gas with N{ = 10, 20, and 30 atoms 
(from bottom to top), trapped in a harmonic potential of 
strength V^/t = 6 X 10 -3 and in a lattice with N s = 100 sites. 
In both panels results of the BALDA/FN scheme at u = +6 
are compared with QMC data. 



suits in the metal-insulator phase-separated regime. At 
u = +2 BALDA shows first indications of a locally in- 
compressible region (a plateau in the density profile at 
n = 1), whereas QMC data still predict the system 
to be fully metallic (compressible). At larger u, the 
plateau also develops in the QMC calculations, but, as 
explained in Sec. lIII A2I BALDA now ceases to converge. 
TLDA calculations, on the other hand, are still possible. 
For fully developed phase separation (cases of u = +6 
and u = -)-8 in Fig. [TTJ TLDA and QMC agree very 
well. At it = +3, the plateau at n = 1 is more pro- 
nounced in TLDA than in QMC, where the density pro- 
file just begins to flatten. In this particular case we 
have also performed a TLDA calculation using the an- 
alytical approach of Ref. 0- Corresponding data are 
labelled by TLDA/LSOC in Fig. El For low densities 
(at the edges of the trap) fully numerical TLDA data 
agree better with QMC, but in the center of the trap 
and at the plateaus TLDA/LSOC and QMC data agree 
better. For the GS energy at this particular value of u, 
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FIG. 8: Top panel: an example of asymmetric trapping 
potential (in units of t) as a function of z; for the case 
V 2 /t = 1.6 x 1CT 2 , V 3 /t = 1.6 x 10 -4 , and V 4 /t = 1.92 x 10 -5 . 
Bottom panel: a double- well potential with V^jt = —4 x 10~ 3 
and Vi/t = 3 x 10" 6 . 



TLDA gives S, 



TLDA 

GS 



16.57* while TLDA/LSOC gives 



^tlda/lsoc _ 25.95t. These numbers have to be com- 



pared with the QMC value £ G 



(16.64 ±0.09)*. 



Using the TF-like TLDA equation l|14fl a simple expla- 
nation for the formation of incompressible regions inside 
the trap in coexistence with metallic regions can be given. 
When the local density reaches unity, i.e. the value asso- 
ciated with the xc Mott gap, the left-hand side of Eq. i|14|) 
takes up the discontinuity A xc (u). This implies that the 
density resists crossing unity and develops instead an in- 
compressible region where the constant value n GS {zi) = 1 
is maintained up to a width W such that the difference 
in the classical potential Un GS (zi)/2 + V cx t(zi), evalu- 
ated at the end points of the incompressible region, ex- 
actly compensates for A xc (u). This criterion allows one 
to find a priori those regions of the trap where the local 
Mott-insulating incompressible phases are energetically 
favorable over metallic phases^. 



V. CONCLUSIONS 

In this work we have shown how a detailed picture of 
ground-state properties of strongly interacting ID ultra- 
cold Fermi gases emerges through a novel DFT scheme 



FIG. 9: Top panel: Site occupation n GS (zi) as a function 
of Zi for a repulsive Fermi gas with Nt — 14 atoms and 
u = +8, trapped in the potential of the top panel of Fig. |S] 
and in a lattice with N s = 50 sites. Results of the TLDA and 
BALDA/FN schemes are compared with QMC data. Bottom 
panel: Site occupation n G s(zi) as a function of Zi for a repul- 
sive Fermi gas with Nf = 20,40, and 60 atoms and u — +2, 
trapped in the potential of the bottom panel of Fig. |S] and 
in a lattice with iV s = 100 sites. Results of the BALDA/FN 
scheme only are presented. 
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FIG. 10: Site occupation n G s(zi) as a function of za for an 
attractive Fermi gas with Nf = 20, 40, and 60 atoms and 
u = —1, trapped in the potential of the bottom panel of 
Fig. |U and in a lattice with iV s = 100 sites. Results of the 
BALDA/FN scheme only are presented. 
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FIG. 11: Site occupation n GS (zi) as a function of Zi for a system of Nt = 70 fermions with u = +2, +4, +6, and +8 in a lattice 
with 7V S = 100 sites, confined by a harmonic potential of strength Vijt = 2.5 x 10~ 3 . BALDA/FN and TLDA results are 
compared with QMC data. 



using as reference fluid a many-body interacting ID 
model which is exactly solvable, the ID homogeneous 
Hubbard model. Needless to say, this basic idea is ap- 
plicable to all inhomogeneous ID systems for which a 
properly chosen, underlying homogeneous reference liq- 
uid can be described by an exactly solvable model, only 
the ground-state energy of the liquid being the key in- 

p^jj .22.23.24.29.46.62 

Our main conclusions are summarized as follows: 

(i) Bethe-Ansatz-based density-functional calculations 
agree quantitatively with independent quantum Monte 
Carlo data. Typical differences between BALDA/FN 
and QMC results for the density profiles in the metallic 
regime are around 1%. In the worst case encountered in 
all the calculations performed in this work the difference 
amounts to 7%. These errors are not easily reduced, but 
are sufficiently small to allow a detailed microscopic de- 
scription of the energetics and the atom-density profiles 
of confined fermions, at a computational cost reduced by 
orders-of-magnitude compared to QMC. Large systems 
of thousands of sites and complex systems, e.g. with 



reduced spatial symmetries, are easily handled by the 
BALDA-DFT scheme. 

(ii) A fully numerical formulation of the BALDA re- 
quires solution of the Lieb-Wu integral equations instead 
of a paramctrization of the energy of the reference fluid. 
The added numerical effort is compensated by a rela- 
tively small gain in accuracy, as judged by comparison of 
BALDA/LSOC and BALDA/FN density profiles to the 
QMC profiles. 

(iii) A total-energy LDA can be set up in the same way 
and provides useful numerical results in situations where 
Kohn-Sham calculations using BALDA for the correla- 
tion energy fail to converge. These situations arise in the 
presence of a local Mott metal-insulator transition, lead- 
ing to phase separation characterized by flat (incompress- 
ible) and metallic (compressible) regions in the density 
profiles. 

(iv) Overall, all four density-functional schemes 
(Kohn-Sham BALDA/FN and BALDA/LSOC where 
they converge, TLDA, and TLDA/LSOC) reveal the 
same physics also seen in the QMC data. DFT schemes, 
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however, tend to predict phase separation (plateau for- 
mation) at slightly lower values of the interaction U. 

(v) Density profiles corresponding to different interac- 
tion strengths, filling factors in the metallic regime, and 
curvatures of the confining potential are quite similar. 
The local compressibility, on the other hand, depends 
sensitively on the details of the system, and can be used 
to discriminate between different choices of system pa- 
rameters (see Fig- EJ- 

(vi) Different external potentials, such as asymmet- 
ric trappings and double- well structures, are easily han- 
dled by the same techniques. Specifically for the double 
well, a signature of tunneling between the two wells is a 
tilted density profile, which shows that atoms can tun- 
nel from one well to the other even for very low filling 
factors. For attractive interactions, additionally density 
waves form separately in each weU^LJS, until the den- 
sity becomes so high that oscillations arising from both 
wells start to merge and a joint pattern develops. This 
type of double- well structure has, to our knowledge, not 
yet been produced optically in ID ultracold atom sys- 
tems, but is readily created in higher-dimensional traps 
and semiconductor heterostructures, to which many of 
our conclusions also apply. 
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APPENDIX A: SOME DETAILS OF SOFT 

The aim of this Appendix is to present a summary of 
the two key results of SOFT—: (i) the Hohenberg-Kohn 
theorem and (ii) the Kohn-Sham mapping to an auxiliary 
noninteracting system. 

The basic variable of SOFT is the site-occupation 
n(zi) = (fy\n(zi)\ty), where is a generic many-body 
state. As in standard DFT, the central result of SOFT is 
the Hohenberg-Kohn (HK) theorem, which can be sum- 
marized in three key statements: (a) the GS expecta- 
tion value of any observable O is a unique functional 
O = (GS|0|GS) = 0[n GS ] of the GS site-occupation 
n GS (zi); (b) the GS site-occupation minimizes the total- 
energy functional £[n]; and (c) £[n] can be written as 

£[n] = ^hkN + V cxt (zi)n(zi) , (Al) 



where JhkM = (vt'|T + 7ii n t|^') is a universal functional 
of the site occupation, in the sense that it does not de- 
pend on the external potential. 

Part (b) of the HK theorem suggests that if the ex- 
act analytical expression of JhrH was known, the GS 
energy and the GS site occupation could be found by 
solving the Euler-Lagrange equation 

^ff^ + = constant , (A2) 

on(Zi) 

the constant having the meaning of a Lagrange multiplier 
to enforce particle-number conservation. 

The Kohn-Sham mapping, again in analogy with stan- 
dard DFT, provides an essential simplification. One con- 
siders a noninteracting auxiliary system described by the 
Hamiltonian 

Hs = -£X>J [4(^)Cafe-)+H.C.] 

+ ^v KS (z t )h(z t ) . (A3) 

i 

The central assertion used in establishing the mapping 
is that for any interacting system there exists a local 
single-particle potential v KS (zi) such that the exact GS 
site occupation n GS (zi) of the interacting system equals 
the GS site occupation of the auxiliary problem n GS (zi) — 

ric} (zi) (noninteracting v- representability) . According to 
part (c) of the HK theorem there then exists a unique en- 
ergy functional £ s [n] = T s [n] v KS (zi)n(zi) , for which 
the variational equation <5£ s [n] = yields the exact GS 
site occupation n G s(zi) corresponding to H s . %[n] de- 
notes the universal kinetic energy functional of noninter- 
acting pseudospin-1/2 fcrmions. 

Suppose that the ground state of H s is nondegenerate. 
The GS site occupation n G s(zi) (and thus, by assump- 
tion, n as (zi)) possesses a unique representation 

n GS (zi) = \¥>a(zi)\ 2 ( A4 ) 

a,occ. 

in terms of the lowest N{ single-particle orbitals obtained 
from the KS-Schrodinger equation 

^ [-ti,j + v KS (zi)dij] (Pa(zj) = e a ip a (zi) . (A5) 

3 

Once the existence of a potential v KS (zi) generating 
n-Gs(zi) via Eqs. (|A4|> and (|A5|) is assumed, uniqueness 
of v KS (zi) follows from the HK theorem. Thus the single- 
particle orbitals Lp a (zi) — y Q [n GS ](zi) are unique func- 
tional of n GS {zi), and the noninteracting kinetic energy 

%[n GS ] = -^2^2^,3 [<Pa{Zi)<Pa{Zj) + C.C.] (A6) 
i,j a 

is a unique functional of n GS (zi) as well. 
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It is convenient at this point to write the total energy 
functional £ G s[ n Gs] m Eq. (|A1|) by adding and subtract- 
ing 7^[n GS ] and a Hartree term £h = Uj^i n Gs( z i)/4, i.e. 



£ G s[™G S ] = %[n GS ] + £H>Gs] + ^Cc[«Gs] 

+ ^2 v o X t(z l )n GS (z i ) , 



(A7) 



where the exchange-correlation functional is formally de- 
fined as £ xc [n GS ] = J^hk^gs] ~ %[n G s] ~ £h[«gs]- Tnc 
HK variational principle ensures that £[n] is stationary 
for small variations Sn(zi) around n GS (z.i): 



£[n GS + 8n] - £ GS [n GS ] = 6% + ^8n{zi) 



V cxt {z t ) + —n GS (Zi) + v xc [n GS ](zi) 



= 0. 



(A8) 



r 



where v xc [n GS ] (zi) denotes the exchange-correlation po- 
tential, 



Vx C [n as ](zi) 



S£ x 



Sn(zi) 



(A9) 



"Gs(zi) 



Using 5% = — '^2 i v KS (zi)Sn(zi), we find that the Kohn- 
Sham potential is given by 

Vns{zi) = vu[n GS ](zi) + v xc [n GS ](zi) + V cxt (zi) , (A10) 

where Va[n as ](zi) = Un GS (z l )/2. 



APPENDIX B: QUANTUM MONTE CARLO 
CALCULATION OF THE GROUND-STATE 
ENERGY 

Within our zero-temperature QMC approach, we use 
the Trotter decomposition 61 in applying the projector op- 
erator to a trial wave function \^t), 

r ~ t«/At 

exp(-6H) |* T > = exp(-ArW) |* T ) ■ (Bl) 



For small values of At one can then split the exponential 
of a sum of noncommuting operators as^i 



exp(-ArW) = exp(-Arr)exp(-ArWint) 

x exp(-Ar7T: cxt ) + 0((At) 2 ) . (B2) 



This is essential for the implementation of the determi- 
nantal QMC algorithm 56 i 58 i 59 i 60 . 

The error introduced by the Trotter decomposition in 
the calculation of the ground-state energy can be shown 
to be of order (At) 2 —. Hence, for our comparison of the 
QMC energies with the ones of the BALDA/FN scheme 
we have made an extrapolation At — > 0. In Fig. 0] 
we show the QMC ground-state energy as a function of 
(At) 2 and the linear fit £§ s mc (At) = £« s mc +A(At) 2 from 
which we get the QMC energies presented in Table [I] 



Electronic address: m.polini@sns.it 

R. Saito, G. Dresselhaus, and M.S. Dresselhaus, Physical 
Properties of Carbon Nanotubes (Imperial College Press, 
London, 1998). 

S.J. Tans, A.R.M. Verschueren, and C. Dekker, Nature 
393, 49 (1998); R. Martel, T. Schmidt, H.R. Shea, T. Her- 
tel, and Ph. Avouris, Appl. Phys. Lett. 73, 2447 (1998); 
Y. Huang, X. Duan, Y. Cui, L.J. Lauhon, K.-H. Kim, 
and CM. Lieber, Science 294, 1313 (2001); B.M. Kim, 
T. Brintlinger, E. Cobas, M.S. Fuhrer, H. Zheng, Z. Yu, 
R. Droopad, J. Ramdani, and K. Eisenbeiser , Appl. Phys. 
Lett. 84, 1946 (2004). 

A. Nitzan and M.A. Ratner, Science 300, 1384 (2003). 
M. Greiner, I. Bloch, O. Mandel, T.W. Hansen, and T. 
Esslinger, Phys. Rev. Lett. 87, 160405 (2001); H. Moritz, 
T. Stoferle, M. Kohl, and T. Esslinger, ibid. 91, 250402 
(2003); T. Stoferle, H. Moritz, C. Schori, M. Kohl, and 



T. Esslinger, ibid. 92, 130403 (2004); BL. Tolra, K.M. 
O'Hara, J.H. Huckans, W.D. Phillips, S.L. Rolston, and 
J.V. Porto, ibid. 92, 190401 (2004). 

B. Paredes, A. Widera, V. Murg, O. Mandel, S. Foiling, 
I. Cirac, G.V. Shlyapnikov, T.W. Hansen, and I. Bloch, 
Nature 429, 277 (2004). 

T. Kinoshita, T. Wenger, and D.S. Weiss, Science 305, 
1125 (2004). 

H. Moritz, T. Stoferle, K. Giinter, M. Kohl, and T. 

Esslinger, Phys. Rev. Lett. 94, 210401 (2005). 

A.H. MacDonald, Phys. Rev. Lett. 64, 220 (1990); X.G. 

Wen, Phys. Rev. B 41, 12838 (1990) and Phys. Rev. Lett. 

64, 2206 (1990); X.G. Wen, Phys. Rev. B 44, 5708 (1991) 

and Int. J. Mod. Phys. B 6, 1711 (1992); for a recent review 

see A.M. Chang, Rev. Mod. Phys. 75, 1449 (2003). 

S. Roddaro, V. Pellegrini, F. Beltram, G. Biasiol, L. Sorba, 

R. Raimondi, and G. Vignale, Phys. Rev. Lett. 90, 046805 



13 



(2003) ; Y.C. Chung, M. Heiblum, and V. Umansky, ibid. 
91, 216804 (2003); S. Roddaro, V. Pellegrini, F. Beltram, 
G. Biasiol, and L. Sorba, ibid. 93, 046801 (2004); S. Rod- 
daro, V. Pellegrini, F. Beltram, G. Biasiol, L. Sorba, R. 
D'Agosta, R. Raimondi, and G. Vignale, Physica E 22, 185 

(2004) ; S. Roddaro, V. Pellegrini, F. Beltram, L.N. Pfeif- 
fer, and K.W. West, Phys. Rev. Lett. 95, 156804 (2005). 

10 R. D'Agosta, R. Raimondi, and G. Vignale, Phys. Rev. B 
68, 035314 (2003); E. Papa and A.H. MacDonald, Phys. 
Rev. Lett. 93, 126801 (2004); R. D'Agosta, G. Vignale, 
and R. Raimondi, ibid. 94, 086801 (2005); E. Papa and 
A.H. MacDonald, Phys. Rev. B 72, 045324 (2005). 

11 F.D.M. Haldane, Phys. Rev. Lett. 47, 1840 (1981). 

12 A. Minguzzi, S. Succi, F. Toschi, M.P. Tosi, and P. Vignolo, 
Phys. Rep. 395, 223 (2004). 

13 W. Kohn, Rev. Mod. Phys. 71, 1253 (1999). 

14 R.M. Dreizler and E.K.U. Gross, Density Functional The- 
ory (Springer, Berlin, 1990). 

15 Density Functionals: Theory and Applications, edited by 
D. Joulbert, Springer Lecture Notes in Physics Vol. 500 
(Springer, Berlin, 1998). 

16 G.F. Giuliani and G. Vignale, Quantum Theory of the 
Electron Liquid (Cambridge University Press, Cambridge, 
2005). 

17 P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964). 

18 W. Kohn and L.J. Sham, Phys. Rev. 140, A1133 (1965). 

19 D. Pines and P. Nozieres, The Theory of Quantum Liquids 
(Benjamin, New York, 1966). 

20 D.M. Ceperley and B.J. Alder, Phys. Rev. Lett. 45, 566 
(1980); B. Tanatar and D.M. Ceperley, Phys. Rev. B 39, 
5005 (1989); F. Rapisarda and G. Senatore, Austr. J. Phys. 
49, 161 (1996); F.H. Zong, C. Lin, and D.M. Ceperley, 
Phys. Rev. E 66, 036703 (2002); C. Attaccalite, S. Moroni, 
P. Gori-Giorgi, and G.B. Bachelet, Phys. Rev. Lett. 88, 
256601 (2002). 

21 F.D.M. Haldane, J. Phys. C 14, 2585 (1981); J. Voit, Rep. 
Prog. Phys. 58, 977 (1995); H.J. Schulz, G. Cuniberti, 
and P. Pieri, in Field Theories for Low-Dimensional Con- 
densed Matter Systems, edited by G. Morandi, P. Sodano, 
A. Tagliacozzo, and V. Tognetti (Springer, Berlin, 2000). 

22 O. Gunnarsson and K. Schonhammer, Phys. Rev. Lett. 
56, 1968 (1986); K. Schonhammer and O. Gunnarsson, J. 
Phys. C 20, 3675 (1987) and Phys. Rev. B 37, 3128 (1988); 
K. Schonhammer, O. Gunnarsson, and R.M. Noack, Phys. 
Rev. B 52, 2504 (1995). 

23 N.A. Lima, M.F. Silva, L.N. Oliveira, and K. Capelle, 
Phys. Rev. Lett. 90, 146402 (2003); K. Capelle, N.A. Lima, 
M.F. Silva, and L.N. Oliveira, in The fundamentals of elec- 
tron density, density matrix and density functional theory 
in atoms, molecules and solids, edited by N.I. Gidopou- 
los and S. Wilson, Kluwer series "Progress in Theoretical 
Chemistry and Physics" (Kluwer, Dordrecht, 2003). 

24 R.J. Magyar and K. Burke, Phys. Rev. A 70, 032508 
(2004). 

25 W. Kohn and A.E. Mattsson, Phys. Rev. Lett. 81, 3487 
(1998). 

26 L.N. Oliveira, E.K.U. Gross, and W. Kohn, Phys. Rev. 
Lett. 60, 2430 (1988). 

27 A. Floris, G. Profeta, N.N. Lathiotakis, M. Liiders, M.A.L. 
Marques, C. Franchini, E.K.U. Gross, A. Continenza, and 
S. Massidda, Phys. Rev. Lett. 94, 037004 (2005). 

28 A. Griffin, Can. J. Phys. 73, 755 (1995). 

29 V.L. Lfbero and K. Capelle, Phys. Rev. B 68, 024423 
(2003); PE.G. Assis, V.L. Lfbero, and K. Capelle, ibid. 



71, 052402 (2005). 

30 S.Y. Savrasov and G. Kotliar, Phys. Rev. B 69, 245101 
(2004). 

31 J.I. Cirac and P. Zoller, Science 301, 176 (2003). 

32 M.B. Dahan, E. Peik, J. Reichel, Y. Castin, and C. Sa- 
lomon, Phys. Rev. Lett. 76, 4508 (1996); E. Peik, M.B. 
Dahan, I. Bouchoule, Y. Castin, and C. Salomon, Phys. 
Rev. A 55, 2989 (1997); B.P Anderson and M.A. Kase- 
vich, Science 282, 1686 (1998); F.S. Cataliotti, S. Burger, 
C. Fort, P. Maddaloni, F. Minardi, A. Trombettoni, A. 
Smerzi, and M. Inguscio, ibid. 293, 843 (2001); S. Burger, 
F.S. Cataliotti, C. Fort, F. Minardi, M. Inguscio, M.L. 
Chiofalo, and M.P. Tosi, Phys. Rev. Lett. 86, 4447 (2001); 
O. Morsch, J.H. Miiller, M. Cristiani, D. Ciampini, and E. 
Arimondo , ibid. 87, 140402 (2001); G. Roati, E. de Mi- 
randes, F. Ferlaino, H. Ott, G. Modugno, and M. Inguscio, 
ibid. 92, 230402 (2004). 

33 D. Jaksch, C. Bruder, J.I. Cirac, C.W. Gardiner, and P. 
Zoller, Phys. Rev. Lett. 81, 3108 (1998). 

34 M. Greiner, O. Mandel, T. Esslinger, T.W. Hansen, and I. 
Bloch, Nature 415, 39 (2002). 

35 M. Girardeau, J. Math. Phys. 1, 516 (1960); E.H. Lieb and 
W. Liniger, Phys. Rev. 130, 1605 (1963). 

36 A. Luther and V.J. Emery, Phys. Rev. Lett. 33, 589 (1974); 
V.J. Emery, in Highly Conducting One-Dimensional 
Solids, edited by J.T. Devreese, R.P. Evrard, and V.E. 
van Doren (Plenum, New York, 1979); EW. Carlson, D. 
Orgad, S.A. Kivelson, and V.J. Emery, Phys. Rev. B 62, 
3422 (2000). 

37 M. Rigol, A. Muramatsu, G.G. Batrouni, and R.T. Scalet- 
tar, Phys. Rev. Lett. 91, 130403 (2003). 

38 M. Rigol and A. Muramatsu, Phys. Rev. A 69, 053612 
(2004) and Opt. Commun. 243, 33 (2004). 

39 M. Machida, S. Yamada, Y. Ohashi, and H. Matsumoto, 
Phys. Rev. Lett. 93, 200402 (2004); M. Rigol, S.R. Man- 
mana, A. Muramatsu, R.T. Scalettar, R.R.P. Singh, and S. 
Wessel, ibid. 95, 218901 (2005); M. Machida, S. Yamada, 
Y. Ohashi, and H. Matsumoto, ibid. 95, 218902 (2005). 

40 X.-J. Liu, P.D. Drummond, and H. Hu, Phys. Rev. Lett. 
94, 136406 (2005). 

41 Gao Xianlong, M. Polini, M.P. Tosi, V.L. Campo, Jr., and 
K. Capelle, cond-mat/0506570 

42 V.L. Campo, Jr. and K. Capelle, Phys. Rev. A , accepted 
and |cond -mat/0508095 

43 W. Hof'stetter, J.I. Cirac, P. Zoller, E. Dernier, and M.D. 
Lukin, Phys. Rev. Lett. 89, 220407 (2002). 

44 See e.g. S. Melinte, M. Berciu, C. Zhou, E. Tutuc, S.J. 
Papadakis, C. Harrison, E.P. De Poortere, M. Wu, P.M. 
Chaikin, M. Shayegan, R.N. Bhatt, and R.A. Register, 
Phys. Rev. Lett. 92, 036802 (2004). 

45 M.T. Bjork, C. Thelander, A.E. Hansen, L.E. Jensen, 
M.W. Larsson, L.R. Wallenberg, and L. Samuelson, Nano 
Lett. 4, 1621 (2004); L. Samuelson, M.T. Bjork, K. Dep- 
pert, M. Larsson, B.J. Ohlsson, N. Panev, A.I. Persson, N. 
Skold, C. Thelander, and L.R. Wallenberg, Physica E 21, 
560 (2004). 

46 M.R. Bakhtiari, M. Polini, Gao Xianlong, M.P. Tosi, V.L. 
Campo, Jr., and K. Capelle, in preparation. 

47 E.H. Lieb and F.Y. Wu, Phys. Rev. Lett. 20, 1445 (1968); 
see also H. Shiba, Phys. Rev. B 6, 930 (1972). 

48 H.J. Schulz, Phys. Rev. Lett. 64, 2831 (1990). 

49 For recent work on 2D electronic quantum dots see e.g. 
M. Gattobigio, P. Capuzzi, M. Polini, R. Asgari, and M.P. 
Tosi, Phys. Rev. B 72, 045306 (2005). 



14 



See e.g. A.N. Kocharian, C. Yang, and Y.L. Chiang, Phys. 
Rev. B 59, 7458 (1999). 

51 J.P. Perdew, R.G. Parr, M. Levy, and J.L. Balduz, Jr., 
Phys. Rev. Lett. 49, 1691 (1982); J.P. Perdew and M. Levy, 
ibid. 51, 1884 (1983); L.J. Sham and M. Schliiter, ibid. 51, 
1888 (1983); W. Kohn, Phys. Rev. B 33, 4331(R) (1986). 

52 N.A. Lima, L.N. Oliveira, and K. Capelle, Europhys. Lett. 
60, 601 (2002). 

53 M. Ferconi, M.R. Geller, and G. Vignale, Phys. Rev. B 52, 
16357 (1995). 

54 G. Sugiyama and S. Koonin, Ann. Phys. (N.Y.) 168, 1 
(1986). 

55 S. Sorella, E. Tosatti, S. Baroni, R. Car, and M. Parrinello, 
Int. J. Mod. Phys. B 1, 993 (1988). 

56 D.J. Scalapino and R.L. Sugar, Phys. Rev. Lett. 46, 519 
(1981); R. Blankenbecler, D.J. Scalapino, and R.L. Sugar, 
Phys. Rev. D 24, 2278 (1981); D.J. Scalapino and R.L. 
Sugar, Phys. Rev. B 24, 4295 (1981). 



J.E. Hirsch, Phys. Rev. B 28, 4059 (1983). 
EY. Loh and J.E. Gubernatis, in Modern Problems in 
Condensed Matter Sciences, edited by W. Hanke and 
Y.V. Kopaev (North-Holland, Amsterdam, 1992), Vol. 32, 
pag. 177. 

A. Muramatsu, in Quantum Monte Carlo Methods in 
Physics and Chemistry, edited by M.P. Nightingale and 
C.J. Umrigar (NATO Science Series, Kluwer, Dordrecht, 
1999), pag. 343. 

F.F. Assaad, in Quantum Simulations of Complex Many- 
Body Systems: From Theory to Algorithms, edited by J. 
Grotendorst, D. Marx, and A. Muramatsu (John von Neu- 
mann Institute for Computing (NIC) Series, Vol. 10, FZ- 
Jiilich, 2002), pag. 99. 
R.M. Fye, Phys. Rev. B 33, 6271 (1986). 
Gao Xianlong, M. Polini, R. Asgari, and M.P. Tosi, 
cond-mat/0511158 



